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Abstract 

We discuss how the vacuum model of Celenza and Shakin with a squeezed 
gluon condensate can explain the existence of an infrared singular gluon prop- 
agator frequently used in calculations within the global colour model. In par- 
ticular, it reproduces a recently proposed QCD-motivated model where low 
energy chiral parameters were computed as a function of a dynamically gen- 
erated gluon mass. We show how the strength of the confining interaction of 
this gluon propagator and the value of the physical gluon condensate may be 
connected. 
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The Global Colour Model (CGM) of Quantum Chromodynamics (QCD), whose main 
aspects have been reviewed in the last years |I[ - ||, is a quark-gluon quantum field theory 
that very successfully models QCD for low energy hadronic processes. In this approach 
an effective gluon correlator models the interaction between quark currents, and quark and 
gluon confinement may appear via the criterion of absence of real q 2 poles for the propaga- 
tors @||]. This is, for instance, the case of an effective gluon propagator with an infrared 
singularity like a delta function S(k) at low energy ||. There are many recent calculations 
exemplifying the remarkable success of this procedure ||[7j • It relates the hadronic properties 
to the Schwinger functions of quarks and gluons, therefore, when comparing the theoretical 
calculations to some low energy data, as pseudoscalars masses and decay constants or other 
chiral parameters, we are learning how is the infrared behavior of the quark and gluon prop- 
agators. As the time goes on this semi phenomenological tool may reveal to be even more 
successful than the relativistic quark model or the bag model. However, the question for the 
mechanism which leads to the infrared singularities present in the gluon propagator in such 
calculations remains open. 

The infrared enhancement of the gluon propagator due to the nonAbelian character of 
the theory and in particular due to the gluon-gluon self coupling, in principle could be 
understood in a rigorous study of the QCD vacuum. Unfortunately the simple perturba- 
tive vacuum is unstable ||, and there is no stable (gauge invariant) coherent vacuum in 
Minkowski space . On the other hand in the context of the construction of a gauge invari- 
ant, stable QCD vacuum in Minkowski space, the squeezed condensate of gluons has become 
a topic of interest to uncover the underlying dynamics of the theory [IIJ - |12 |. Within this 



class of vacuum models, Pavel et al. |L3| recently proposed a phenomenological vacuum based 



on Abelian QCD which leads exactly to the infrared singularity in the gluon propagator as 
considered in Ref. ||. On the other hand calculations within the GCM suggest that the 
best fit of chiral symmetry breaking parameters are obtained with a propagator containing 
a delta function plus a propagator that behaves as 1/k 2 in the ultraviolet (consistent with 
QCD) but damped at k 2 = Unfortunately, this is not the case of Ref. |L3|] where a 



simple 1/k 2 is obtained together with the delta function. 

Even if a model of the squeezed QCD vacuum is not determined from first principles, its 
properties may be very representative if they lead to a consistent phenomenology, indicating 
the path to the actual vacuum. In this note we show that the model of the QCD vacuum 
proposed by Celenza and Shakin |TD| , |TTf reproduces completely the gluon propagator of 



Ref. 0, i.e. it gives a singular part a la Munczek and Nemirovsky [f| as well as a piece 



containing a dynamically generated mass fll5 , |l6 |. The effective dynamical gluon mass is the 



unique parameter in the model, and the gluon propagator obtained according to Ref. fllOlJll 
is totally compatible with the idea of the GCM. We will briefly discuss the CGM and 
justify a gluon correlator which reproduces many aspects of the chiral symmetry breaking 
phenomenology. Then, we show that this gluon correlator naturally appears in the squeezed 
vacuum model of Ref. |T0|JTT[] , verifying that the parameters involved in the model (the 



strength of the confining interaction, the dynamical gluon mass and the value of the physical 
gluon condensate) are consistent with the ones in the literature. The coincidence between 
the models may indeed suggest an interesting role for the Celenza and Shakin vacuum model. 

The action of the GCM can be obtained from the QCD generating functional through 
the standard method presented in Ref. |JXj] : 

Z[fj, 7]} = N J DqDqDAexp (-S[q, q,A} + J d 4 x(fjq + qr})j , (1) 

where 

S[q, q,A} = J d'x (q{$) + m- ig^ A a )q + \(P^G%\ , (2) 

and £r* = d^A" — d u A® + gf abc A b ^A c u . In the above equation we have not written the gauge 
fixing term, the ghost field term and its integration measure. Introducing a source term for 
the gauge field and writing 

exp(W[J]) = J DAexp (- J d'x^G^G^ - JJA")) , (3) 

the generating functional becomes 



Z[t}, rj] = N J DqDqexp(—q(0 + m)q + r/q + qrj)exp \W 



W~2 If i*Q 



(4) 



where the spacetime integration is implied. The functional W[J] is the generator of con- 
nected gluon n-point functions without quark-loop contributions. It may be written as 

W[J] = \J d A xd 4 yJ a »(x)g 2 D, v (x - y)J a u (y) + W R [J}. (5) 

The main characteristic of the GCM is the fact that we neglect the higher n-point functions 
contained in Eq.([5]) and expressed by H^[J]. The effect of this approximation can only be 
measured in the model building, but it is expected that the phenomenological propagator 
g 2 D llv (x — y) retain most of the information about the non-Abelian character of QCD. With 
this approximation the generating functional can be factorized as 



S X a 



> 'GCM 



hv\- ( 6 ) 



Z[f},rj] = exp \ W R 
Z GC m is giving by 

Z G cm[Vi r}] = N f DqDqexp (-S GC m% q]+fjq + qrj) , (7) 

with 

S GCM % q] = J d A xq{x){$) + m)q(x) + ±J d A xd 4 yJ^x)g 2 D^(x - y)J a v {y). (8) 

The action Sgcm[q, q] together with the generating functional Z GC m[Vi v] defines the GCM. 
The idea of the model is that the nonperturbative behavior that could be missed in the 
truncation performed above can be mostly represented by an effective model in the infrared 
of the gluon propagator in Eq.@. 

There are several discussions in the literature about the nonperturbative behavior of 
the gluon propagator [fJ,|Tj],[II|] , as well as ansatze motivated by an impressive fitting of the 
low energy QCD phenomenology. One such case is the propagator proposed by Frank and 
Roberts which yields the expected QCD behavior in the ultraviolet and presents an 
integrable singularity at the origin. This is accomplished in Landau gauge by the following 
form (in the sequence all the momenta are in Euclidean space) 



k 2 [l + U(k 2 )] 



(9) 



where, 



A(A; 2 ) = 



k 2 [l + U(k 2 )] 



4n 2 d 



ATT 2 m 2 t 5 \k) + 



; (-fc 2 /[4m?D 



(10) 



with d = 12/(33 — 2nf), and rif = 3 (considering only three quark flavors). The mass scale 
rrit determined in Ref. [0] was interpreted as marking the transition between the perturbative 
and nonperturbative domains. We recently modified the above gluon propagator in order to 
corroborate with the suggestion of several works [|l^,|l(| that the gluon can develop a dynam- 
ical mass and established that the unique parameter, m t , could conveniently be substituted 
by the scale associated to this mass. Then, the Eq.fllOl) should be modified to ]14] 

,2 



A(k 2 



k 2 [i + n(fc 2 )] 



4n 2 d 



^ 2 m 2 g b\k) + 



1 



k 2 + m 2 {k 2 ) 



(11) 



where, 



m 2 (k 2 ' 



m 



9 k 2 + m 2 



(12) 



is the dynamical gluon mass. Eq.(^) interpolates between the large momenta behavior of 



the gluon mass predicted by the operator product expansion [18 

34A% 2 l^G^G 



m 2 (k 2 ' 



(13) 



9(N 2 - 1) k 2 

and its finite constant infrared behavior as discussed in Ref. ||15|| . In Eq.(|T3|) i^-G^G^ 
is the gluon condensate. Eq.(|TT]) has in its first term the singularity of Eq.(^), and the 
dynamical gluon mass gives a natural (and needed) damping of the second term at small 
k 2 . In Ref. we obtained a good fit to some chiral parameters ( as the pion mass and 
the quark condensate) using the value m g ~ 600 MeV. Our next step is to show how the 
model of Ref. [|T0| , [TT]1 can fully describe the behavior of Eq . (]TT|) . 



Colour- singlet coherent states can be formed integrating over the group elements associ- 
ated with the colour group 

|Z„(t) >=N J[dg]U(g)\Z(t) >, (14) 

where N is a numerical factor, g specify an element of the group SU(N) and U(g) are 
rotation operators whose properties are described in Ref. [IT|. Celenza et al. [TTJ were able 
to form a squeezed state |Z G > such that 

< Z | : E a (r, t) ■ E a (r, t) : |Z >= 0, (15) 

and 

y < Z | : B (r, t) • B (r, t) : |Z >^ 0. (16) 

Meaning that the condensate is purely magnetic (B a (r, t)), while the color-electric field 
(E a (r, £)) leads to a vanishing condensate. 

Many of the properties of this model were obtained previously with the assumption that 
the gluon field could be decomposed into a constant condensate field and the quantum 
fluctuations g®(x) around it |jT0fl , 



G>(*) + (17) 

Taking into account Eq. fli"7|) and < gt{x) >= 0, one obtains the decomposition of the 
nonperturbative gluon propagator into two parts, 

g*G*{x - y) =< g 2 QlQ h v > + < g 2 g%x)g h v {y) > . (18) 

Note that the choice of Eq. ([L7]) is a purely phenomenological one, and the correct construc- 
tion of the model goes through the steps of Ref. fll|| . However, in this way it is easier to 
verify some of the consequences of the model. First, it does generate a massive effective 



lagrangian [K| 



1 „„ , s , ml 



C mg (x) = -- G ; u (x)G^(x) + -fgZ(x)gr u (x) + ... (19) 



Secondly, the gluon mass m 2 is equal to (15/32) < g 2 Q 2 >. In Ref. [[Hj the square of this 



condensate has been identified with the phenomenological gluon condensate of Shifman et 
al. |P^| , but we do not pursue this point since the gluon mass may contain a large part 
of dynamics, and we do not see any a priori reason for the identity of the condensates 
(although it looks phenomeno logically plausible). Finally, what more interest us is that in 
Euclidean momentum space the effective nonperturbative gluon propagator corresponding 
to the decomposition of Eq.(|l8|) has the form |f20| 



g 2 A(k 2 ) oc 



1 k 2 + m 2 g 



(20) 



where the S 4 (k) comes from the condensate field and the massive propagator is a consequence 
of both fields (since the condensate field generates the gluon mass). Once we have not 
saturated the condensate resulting from Eq.(|18D with the one of Shifman et al., the value 
of the dynamical gluon mass can be considered a free parameter and not yet related to the 
gluon condensate. Note that we do have some freedom in the definition of the factors of 
Eq.flTTD []14[ , as well as in the definition of the gluon mass in Eq. (PD|) , this is why we may 



affirm that Eq.flTTD an d Eq.(|2"0"D are totally compatible. The important point is that the 
form of the propagator in Eq.fl2"0T) is consistent with the propagator of Ref. and gives 
a meaning for the confining propagators of the GCM. It is interesting that the proposal of 
Ref. was to modify the propagator of Ref. introducing the concept of the dynamical 



gluon mass, and it turned out that this modification just led to the model of Ref. [|10] 



In the Celenza and Shakin model the gluon mass and the strength of the confining 
interaction are related to the constant condensate field, and the square of this one was 



identified with the gluon condensate of Shifman et al. |TJ|. As we said above it is not 
clear to us how much of the dynamics can be built over this condensate, and we can let 
the gluon mass as a free parameter in Eq. (|20|) . However, we would like to show that this 
identity is natural. There are several ways to verify that the gluon mass (or the condensate 
field) are related to the phenomenological gluon condensate. We obviously expect that the 
dynamical gluon mass is connected to the gluon condensate through the operator product 



expansion as described by Eg. (|T3|) . We can also compute the vacuum energy as a function 
of the propagator in Eq. (|20|) , and than equalize this vacuum energy to the trace anomaly 
expression 4^ (G^G^), where (3(g) is the perturbative f3 function of the renormalization 
group equation. This approach was used in the last paper of Ref. [[L5[ and more recently in 



Ref. [IT], where the relation with the gluon condensate was explored to fix the gluon mass. 
As these methods have already been discussed for massive gluons, we propose here other 
procedures where we can show that the gluon propagator of Eq.(^) can indeed be related to 
the gluon condensate as well as to other typical hadronic parameters. We start by showing 
the relation between the confining interaction and the bag constant. 



The bag constant (B) for a given gluon propagator can be obtained through |2T] 



rk 2 dk 2 ( ( A 2 (k 2 )k 2 + B 2 (k 2 ) \ B 2 (k 2 ) \ 

J (2n)*\ \ A 2 (k 2 )k 2 ) A 2 (k 2 )k 2 + B 2 (k 2 ) J ' 1 1 

where A(k 2 ) and B(k 2 ) appear in the inverse of the renormalized quark propagator 

S~ l (k) =ift + E(jfe) = i fA(k 2 ) + B(k 2 ). (22) 

The form of A(k 2 ) and B(k 2 ) is obtained solving the Schwinger- Dyson equations for a given 
gluon propagator, and here enter the information about the GCM propagators. 

The Schwinger-Dyson equations for a gluon propagator equal to Eq.(|TTl) were discussed 



in Ref. [£L4j] . To avoid a numerical calculation that may conceal the simplicity of the problem, 
we will approximate the gluon propagator only by its confining part. The part proportional 
to (k 2 + m 2 )" 1 is responsible for the tail of the propagator, actually, the larger is the gluon 
mass the smaller is the contribution to this tail. Therefore, we will not introduce a large 
error in neglecting the second term of Eq.flTTD or Eq. fl20|) . Proceeding in this way it is easy 
to verify that the propagator A(k 2 ) = lQTT i m 2 5 A (k) imply in the following solution for A(k 2 ) 



A(k 2 ) = 2 , B(k 2 ) = 2Jm 2 - k 2 . (23) 



and B(k 2 ) when k 2 < m 



For k 2 > m 2 we have B(k 2 ) = 0. With the solution of Eq.fl2"3"|) the calculation of B is 
straightforward and gives 



For a gluon mass of approximately 600 MeV |14jl5| we obtain B = (169 Me V) A which, 
considering the approximation performed to obtain Eq.(^), is in good agreement with the 
MIT value of (UQMeV) 4 . 

The bag constant can be related to the string tension for fermions in the fundamental 



representation through 22 



K F = (87ra s C F B) 1/2 , 



(25) 



where Cp is the quadratic Casimir operator for the fundamental representation. On the 
other hand this same string tension has been estimated in the last paper of Ref. |T^] as 



7T" \ ^fiu 



mi 



(26) 



From the identity of Eq. (1251) and Eq. (|2"6|) we obtain 



7T / :j 
9 V27ra« 



GTG 



fJIV 



(27) 



Eq.(p7|) shows the connection between the gluon mass and the gluon condensate. Taking 



1 and l^G^G 



(0.01) GeV [19| we obtain m g w 550 MeV, which is consistent 



with the value we quoted before (see Refs. | 14],|10 |). 

Another way to verify the relation between the gluon mass and the gluon condensate is 
comparing the bag energy density to minus the vacuum expectation value of the trace of 
the QCD energy momentum tensor, although it is known that this comparison does not give 



reasonable numbers |]23| . The energy density in the bag model at the equilibrium radius is 
equal to 



n = 4B, 



(28) 



and from the vacuum expectation value of the trace of the QCD energy momentum tensor 
we obtain 



1 1 n 

= l (e w > = -~( llN - 2n f ) (G^CT) . (29) 

which gives 

^= (11JV - 2 "^ 2 (^ G ^). (30) 

As expected we obtain a gluon mass proportional to the gluon condensate, but about 1.5 
times the desired value (considering rif = 3). 



In conclusion, we discussed how the QCD vacuum model of Celenza and Shakin [11 



explain the presence of a confining gluon propagator frequently used in calculations within 



the global colour model. In particular, it reproduces the propagator of Ref. |14] where the 
dynamically generated gluon mass plays a fundamental role. We considered the strength 
of the confining interaction, equal to the gluon mass, to be a free parameter contrarily to 



Ref. [rU|, based on the fact that the dynamics may modify the gluon propagator. We have 
given arguments showing that the gluon mass is indeed related to the phenomenological 
gluon condensate, as well as it can be related to the string tension and bag constant. The 
coincidence between the propagators used in applications of the global colour model and the 
one resulting from the squeezed vacuum model of Celenza and Shakin might indicate some 
intrinsic property of the actual QCD vacuum. 
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